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The factorization theorem in decays of B( 3 ) mesons to two charmed mesons (both pseudoscalar 
and vector) can be proved in the leading order in mo/mB and Aqcd/?tid expansion. Working 
in the perturbative QCD approach, we find that the factorizable emission diagrams are dominant. 
Most of branching ratios we compute agree with the experimental data well, which means that the 
factorization theorem seems to be reliable in predicting branching ratios for these decays. In the 
decays of a B meson to two vector charmed mesons, the transverse polarization states contribute 
40% — 50% both in the processes with an external W emission and in the pure annihilation decays. 
This is in agreement with the present experimental data. We also calculate the CP asymmetry 
parameters. The results show that the direct CP asymmetries are very small. Thus observation of 
any large direct CP asymmetry will be a signal for new physics. The mixing induced CP asymmetry 
in the neutral modes is large. This is also in agreement with the current experimental measurements. 
They can give a cross check of the sin 2/3 measurement from other channels. 

PACS numbers: 13.25.Hw, 12.38.Bx 



I. INTRODUCTION 

The hadronic decays of B meson are important for particle physics since they provide constraints of the standard 
model Cabibbo-Kobayashi-Maskawa (CKM) matrix, a test of the QCD factorization, information on the decay mech- 
anism, and the final state interaction. The CP asymmetries, in which some of the hadronic uncertainties are canceled 
in their theoretical predictions, play an important role in the investigations of B physics. For the decays with a single 
D meson in the final states, only tree operators contribute, and thus no CP asymmetry appears in the standard model 
However, for decays with double-charm final states, there are penguin operator contributions as well as tree 
operator contributions. Thus the direct CP asymmetry may be present. Recently, the Belle Collaboration reported a 
large direct CP violation in B° — > D + D~ decay Q, while BaBar reproted a small one, with a different sign even Q]. 



What is more, large direct CP asymmetries have not been observed in other B° —> D* + D* decays 



3 ei 



either, which 



have the same flavor structures as B' —> D + D at the quark level. Intrigued 



by these experimental results, many 



investigations on the decays of B to double-charm states have been carried out 

The theoretical study of hadronic B decays has achieved great success in recent years. Among them, the perturbative 
QCD approach (PQCD) is based on fcy factorization Q. By keeping the transverse momentum of quarks, the end 
point singularity in the collinear factorization has been eliminated. Since transverse momentum introduces another 
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energy scale, double logarithm appears in the QCD radiative corrections. The rcnormalization group equation is used 
to resum the double logarithm, which results in the Sudakov factor. This factor effectively suppresses the endpoint 
contribution of the distribution amplitude of mesons in the small transverse momentum region, which makes the 
perturbative calculation reliable. Phenomenologically, the PQCD approach successfully predict the following: (l)the 
direct CP asymmetry in B decays [l(3], (2)the pure annihilation type B decays 11| (3) the strong final state interaction 
phase and color suppressed decay amplitude in the B — > Dir decays [lj. 

In charmless two-body B decays, the final state mesons can be considered as massless therefore both of the final 
state mesons are on the light cone. The collinear factorization can be easily proved in the heavy quark limit. For the 
decays with a single heavy D meson in the final states, one can still prove factorization 12( in the leading order of the 



r = mu/rriB expansion. For the decays with double-charm quarks in the final states, such as B — > J/tpK 1 XcK, it is 
believed that the factorization fails. However the decays with double D mesons in the final states are different. The 
reason is that the expansion parameter mo/ms ~ 0.36 could be considered small, rnj/^/ms ~ 0.6 is not. In other 
words, the J/ip (xc) are soft particles in B decays; while the meson is collinear in the B — > D^D^ decays. 
The momentum of the Dj?} meson in the latter decays is |p| ~ |ms(l — 2r 2 ), which is still nearly half of the B meson 
mass. The decays of B to double-charm states can be investigated in the PQCD approach in the leading order of 
r = mn I m-B & n d Aqce> /™d expansion. All of the annihilation type diagrams contain end-point sigualrity, which are 
quite different from the spectatorlikc diagrams which are dominated by the form factors. It is very difficult to deal 
with in the collinear factorization. The PQCD base on kx factorization is almost the only approach that can give 
quantitative calculations of annihilation type decays. 

This paper is organized as follows. In Sec. UH we list the formalism, including the Hamiltonian, the wave functions 
of the mesons, the factorization formulae of the Feynman diagrams for B — > PP decay mode, and the analytic 
expressions for the decay amplitudes. In Sec. IIII1 the numerical results of the physical obscrvablcs and discussions 
of the results are given. Sec. IIVI is a brief summary. The common PQCD functions, scales, and the factorization 
formulae of the Feynman diagrams for B — > PV, B — > VP, and B — > VV modes are all put into the appendices for 
simplicity. 



II. ANALYTIC EXPRESSIONS 



In hadronic B decays there are several typical energy scales, and expansions with respect to the ratios of the scales 
are usually carried out. The physics with a scale higher than the W boson mass are electroweak interactions, which 
can be calculated pcrturbativcly. The physics between the W boson mass and b quark mass obtain QCD corrections. 
This correction is included in the Wilson coefficients of the four-quark operators in the effective Hamiltonian. The 
physics below the b quark mass is more complicated. We have to utilize the factorization theorem to factorize the 
nonperturbative contributions out, so that the hard part can be calculated pcrturbatively. In the PQCD approach, 
we utilize the kx factorization where the transverse momenta of the quarks in the mesons are kept to eliminate 
the end-point sigularity. Because of the new transverse momentum scale introduction, double logarithms appear 
in the calculation. We resum these logarithms to give a Sudakov factor, which effectively suppresses the end-point 
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region contribution. Thus the end-point sigularity in the usual collinear factorization disappears. This makes the 
perturbative calculation reliable and consistent. For decays with D meson in the final states, another scale mo is 



introduced. The factorization is proved in the leading order of the m d /tub expansion 12| , therefore, as it is done in 
the computation of B — > DM and B — > DAI amplitudes [3] , we will work in the leading order m^/ms expansion. 
For each of the diagrams in the following, we keep the contributions in the leading order of mo /tub- For example, in 
the B meson to two vector mesons decays, the leading order contributions of some transversely polarized amplitudes 
are proportional to r 2 (r = TOo/ms). Then we will keep the r 2 terms in these diagrams. While in other cases, the 
terms of r 2 are neglected because the leading order is lower than 2. Finally the amplitude for B — > M2-M3 (M2 and 
M3 stand for two mesons) decay within PQCD approach is decomposed as 

M = y"rf 4 fc 1 d 4 fc 2 d%$ B (fc 1 ,t)T ff (fc 1 ,fc 2 ,fc3,i)$M 2 (fe,t)$M 3 (fc3,0e 5(fc - t) , (1) 

where fcj (i = 1, 2, 3) are the momenta of the quarks in mesons which are defined explicitly in Eq.([10|). Th is the hard 
part that is perturbatively calculable. 3>b and (i = 2, 3) are the hadronic meson wave functions that are treated 
as nonperturbative inputs. The Sudakov factors e^*"^*) (i = 1 7 2, 3) are from the rcsummation of double logarithms . 



A. Notations and conventions 



The Hamiltonian referred to in this paper is given by \h\ : 

n ^ = %{ E V qb V* D [C 1 { l j)0l^) + 02(^)01^) 

* n—11. r 



q—u,c 

10 . 

-V*V? D [X>M W] \ + H.c, (2) 

i=3 ' 

where V q b{D) and V tb (D) with D = d,s are CKM matrix elements. Functions Oi (i = 1, 10) are local four-quark 
operators : 

• current-current (tree) operators 

0\ = (q a bp)v-A(Dpq a ) v ~A, 0\ = (q a b a )v^A(Dpqp) V - A , (3) 

• QCD penguin operators 

3 = (D a b a ) v -A^2(q' p q'p)v-A, 4 = (D p b a )v-A E(?U/;)v-A, (4) 

5 = (D a b a ) V - A ^2(q'l3q'p)v+A, 6 = {Di 3 b a ) v ^ A ^2{q' a q' p ) v+A , (5) 

q' q> 

• electro- weak penguin operators 

3 3 

7 = -(Dab^v-A^eg'iq'pq'^v+A, 8 = -{Dpb a ) V - A E e q' (vWph+A, (6) 

q' q< 

3 ^ 3 

O9 = -{D a b a ) V -A^2e ql (q' p q'p) V -A, O10 = -{Dpb a )v-A E e 9'(9a^)v-A, (7) 



where a and /3 are color indices and q' are the active quarks at the scale m&, i-e. q' = (u,d, s,c,b). The left-handed 
current is defined as (q' a q'p)v-A = q' a lv{^- - l5)q'p and the right-handed current is {q' a q'p)v+A = q' a lv{l + Ibjq'p- The 



combinations a, of Wilson coefficients are defined as usual 
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01 = C 2 + Ci/3, a 2 = d + C 2 /3, a 3 = C 3 + C 4 /3, a 4 = C 4 + C 3 /3, a 5 = C 5 + Gs/3, 

a 6 = C 6 + C 5 A «r = C 7 + C 8 /3, a 8 = C 8 + C 7 /3, a 9 = C 9 + C 10 /3, a 10 = C w + Cfe/3. (8) 

We work in the light-cone coordinate, in which a vector is defined as ( V °^ 3 ; ^'j^ 1 V 1 ^ 2 ). We use Af 2 to 
denote the charmed meson with a c quark and A/3 to denote the meson with a c quark. In this paper we work in the 
rest frame of B meson and define the direction in which M 2 moves as the positive direction of z-axis. Therefore the 
momenta of Bi s \ meson and two charmed mesons are defined in the light-cone coordinate as 

PB = 1, ± ), p 2 = - rlrl 0_l), p 3 = ^L(rl 1 - r 2 , 0± ), (9) 

where r,; = mi/rrig (i = 2,3) and 0j_ are zero two-component vectors. to 2 and m 3 are the masses of the two 
charmed mesons. One can find that our definitions of the momentums violate the on shell conditions. In the following 
calculations we will keep the contributions of each diagram to the leading power of r^(i = 2,3). One will find that 
all the terms with a power of r, higher than 2 arc dropped. At this accuracy level, the on-shell conditions can be 
satisfied. We use k\, fc 2 , and k 3 to denote the momenta carried by the light quarks in Bt s \ meson and two charmed 
mesons. They arc defined by 

fcl = (0 ,^| Xl ,k 1± ), fc 2 = (^(l-r 3 2 )x 2 ,0,k 2J J, fc 3 = (0,^(l-r|)* 3 ,k 3± ), (10) 

with xi, x 2 and x 3 as the momentum fractions. 



B. Wave functions of B( s i mesons 



The -B(s) meson wave functions are decomposed into the following Lorentz structures: 

— - ttl -*<0|M0)d fl (*)|iJ w (Pi)) 



(2rr)4 



Pa 

Here, 4>B {3) {k\) and <t>B, r) (fci) are the corresponding leading twist distribution amplitudes, and numerically (f>B {s) {ki) 
gives small contributions [16|, so we neglect it. The expression for 3>b (s) becomes 

*B W = -y=(/Pl + m S (a ,)750i3 (s) (fcl). (12) 

For the distribution amplitude in the b-space, we adopt the model function 



4>B (B) {x,b) = N B[s) x 2 (l - x) 2 exp 



1 xm Bi , 

2 ujb 



(13) 
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where b is the conjugate space coordinate of ki^. Nb {s) is the normalization constant, which is determined by the 
normalization condition 

fi 



h 



(14) 



The and B® decays are studied intensively in PQCD approach[9j. With the rich experimental data the tub = 



0.40 GeV is determined for B meson. For B s meson, we will follow the authors in Ref. 
u ba = (0.50 ±0.05) GeV. 



28[ and adopt the value 



C. Wave function of meson 



In the heavy quark limit, the two-particle light-cone distribution amplitudes of £)W/Z)W meson are defined as 17 



(D(P 2 )\q a (z)cp(0)\0) = 



(D*(P 2 )\q a (z)cp(0)\Q) 



1 



dxe 



ixPi -z 



[7s(/*2 + m D )4> D (x,b)] aj3 



c Jo 



dxe 



IXP2 -z 



[fr{p2 + m D >)<j> L D *{x,b)+ MA+m D ')cf>l,(x,b)] (15) 



with 



fT 



as the normalization conditions. In the heavy quark limit we have 



m c + m d 



T m c + m d 



(16) 



(17) 



Thus we will use = fn* in our calculation. The model for the distribution amplitude for D meson that we used 
in this paper is 



4> D {x,b) 



2J2N, 



f D 6x(l -x)[l + C D (1 - 2x)] cxp[^^] 



(18) 



which has been tested in the B D^M and B -> D^M decays [l3|]. The masses of D, ! meson that we use are 



m D = 1.869 GeV, 
m D , = 2.010 GeV, 



1.968 GeV, 



m D ,- = 2.112 GeV. 



(19) 



We use C D = 0.5 ± 0.1, uj = 0.1 GeV for DID meson and C D = 0.4 ± 0.1, u) = 0.2 GeV for D s /D s meson, which are 
by fitting. In the wave function of Df ■, mesons, the 4>j^, and (jr£>, can not be related by the 



equation of motion. We simply follow the authors in Ref. 



0) 

and adopt the same model as that of D meson for them 



0h*(x,b) = 4>l,(x,b) = 



f D .6x(l -x)[l + C D .(1 - 2x)} exp[ - ] 



2J2N, 



(20) 



The mass difference of Z?( s ) and D*^ is very small. In a heavy quark limit, the light meson in D) s ( mesons is not 
sensitive to the spin and color of the heavy c or c quark. Thus the light-cone distribution amplitudes of Dr s \ and 
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. C 

M 2 



(a) (b) ( c ) (d) 



FIG. 1: Emission diagrams. 



DJs should be very similar. In our calculation, we simply take Co* = Cd and lu* = cj. fn = (207 ± 4) MeV 21 1 
and fj} s = (241 ± 3) MeV[21| are adopted and the following relations derived from HQET [22( are used to determine 
Id* 



(») 



Id, 



f Dl - = 



"D7 



for- 



(21) 



The value of fo s above is smaller than the recent experimental data Jd s = (273±10) MeV [18|. Because the amplitude 
in the PQCD approach is factorized as the convolution of the wave functions, Sudakov factors and the hard part, the 
branching ratio is proportional to the /m 2/3 - Thus if the experimental data is adopted, our results for the branching 
ratios will increase by F = ( ^rnp ) 2 for single D s meson in the final state and F 2 for double D s meson final state. 



D. Factorization Formulae for B — > PP mode 



In this subsection, we list all the amplitudes from the Feynman diagrams for (M2Mz\Ci(ii)Oi{ii)\Br s \) up to the 
leading order, with Mi and M3 as two charmed mesons. According to their topological structures, the diagrams that 
contribute to the decays of to two charmed mesons can be divided into two types, the emission diagrams (see 
Fig. [T]) with the light antiquark in Bi s \ meson entering one of the charmed mesons as a spectator and the annihilation 
diagrams (see Figs. [2] and [3]) without any spectator quark. The first two diagrams in Fig.[T]arc factorizable diagrams, 
whose amplitude can be naively factorized as a decay constant of a charmed meson and a form factor like structure. 
The amplitudes arise from all the possible Lorentz structure of the operators for factorizable emission diagrams are 
given as following, where a* denotes the Wilson coefficients and t is the scale. 

• Factorizable emission diagrams for (V-A)(V-A) operator 

,1 ,1/A 

F^ L {ai{t)) = 8irC F fM 3 m 4 B / dxidx 2 \ bidbib 2 db 2 (t)B{xi,bi)4>M 2 (x2) 



E e {tP)a i (t^)h e (x 1 ,x 2 {l - r 2 )MM)St{x2){l + x 2 + r 2 (l - 2x 2 )) 
r 2 {\ + r 2 )E e {t^)a l {t^)h e {x 2 ,x 1 {\~rl)MM)St{xi)\ . (22) 



-Hi 



® 



M 2 
(a) 



® 



Ma 
(b) 



M 
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FIG. 2: Annihilation diagrams without charm quark in the four-quark operator. 
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(a) 



c A/, 



® 



A/-. 



.1/-, 



A/ 3 



M 3 
(c) 



-1/:; 



FIG. 3: Annihilation diagrams with charm quark in the four-quark operator. 



Factorizablc emission diagrams for (S-P)(S+P) operator 



rl ,1/A 

F e p (ai(t)) = l6TrCFfM 3 m% / dx\dx 2 / &i dbib 2 db 2 (j)B{x\, bi)4>M 2 {x 2 ) 

Jo Jo 

E e (tP)a 1 {tP)h e (x 1 ,x 2 (l - r 2 3 ),b ll b 2 )S t (x 2 )(l + 2r 2 + r 2 x 2 ) 



xr 3 



+r 2 E e (t^)a t (t^)h e (x 2 ,x 1 (l - r§), b 2 , b^Sfa) 



The amplitudes for the nonfactorizable emission diagrams in Figfljc) and (d) are given as: 

• Nonfactorizable emission diagrams for (V-A)(V-A) operator are 

l~2 r 1 r 1 / A 

F^ n L (ai(t)) = I6iv\ -C F m B / [dx] \ b 1 db 1 b 3 db 3 (\) B {x 1 M) ( i>M- i {x , 2)^M 3 {x 3 ) 

V o J Jo 

{x 3 - r 2 x 2 )E b {t^)ai{t^)h^{xi,bi) 
(x 2 r 2 - x 2 + x 3 - l)E m (tW)ai{tW)hW(xi, h) 



Nonfactorizable emission diagrams for (V-A)(V+A) operator are 



\l —Cpm 



1 f l/A 

% j [ax\ 
Jo 



xr 3 



[dx] / bidbib 3 db 3 (/)B{xi,bi)(j)M 2 (x 2 )(i)M 3 (x 3 ) 
Jo 

(x 3 + r 2 {x 2 + x 3 ))E en (t^)a i (t^)h'^(x i ,b i ) 



-(r 2 (x 2 -x 3 + 2)-x 3 + 2)E en (t^)a i {t^)h^(x i ,b i ) 



(23) 



(24) 



(25) 



The first two diagrams in Figs. [2] and [3] are the factorizablc diagrams for annihilation diagrams, whose amplitudes 
can be factorized as a £?( s ) meson decay constant and a form factor like structure between two charmed mesons. It 
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should be reminded that, in the decays we considered, the factorizablc annihilation diagrams can be divided into two 
types, depending on whether the quark propagator is a light quark propagator (see the first two diagrams in Fig. [2]) 
or a c-quark propagator (see the first diagrams in Fig. [3]). In calculation we keep the mass of the c-quark while the 
mass of the light quark is neglected and thus these two types of diagrams have different expressions. The amplitudes 
for the factorizablc annihilation diagrams with a light quark propagator (the first two diagrams in Fig. [2]) are given 
as follows: 



• Factorizablc annihilation diagrams for (V-A)(V-A) operator 

»1/A 



b 2 db 2 b 3 db 3 cj>M 2 {x 2 )(j>M 3 (^3) 



Fa L (Mt)) = 87rC F f B m% f dx 2 dx 3 f 

Jo Jo 

x [(2r 2 r 3 (x 2 -2)+x 2 -l) 

xE a (t^) ai (t^)h a (l - (1 - r 2 2 )x 3l 1 - (1 - rl)x 2 MM)St{x 2 ) 
+ (-2r 2 r3(x3-2)-(x3-l)) 

xE a (tW) ai (tW)h a (l - (1 -r 2 )x 2l l - (l-r 2 )x 3 ,b 2 ,b 3 )S t (x 3 )} . (26) 



The two terms of (eij(i)) give destructive contributions. Very little contribution appears when (f)M 2 an d 



KM,, 



are different from each other. Otherwise, F^ L (ai(t)) = 0. 

• Factorizable annihilation diagrams for (V-A)(V+A) operator 

F a LR (a t (t))=F a LL (a,(t))- (27) 

• Factorizable annihilation diagrams for (S-P)(S+P) operator 

,i/A 



Ff P Mt)) = l67rC F f B m% f dx 2 dx 3 [ 

Jo Jo 



b 2 db 2 b 3 db 3 <j) Ma (x 2 ) 4>m 3 (x 3 ) 

x (2r 3 + ra(l - x 2 ))E a (t^) ai {t^)h a {l - (1 - r 2 2 )x 3 , 1 - (1 - r 2 3 )x 2 ,b 3 , b 2 )S t (x 2 ) 

+ (2r 2 + r 3 (l - x 3 ))E a (t^) ai (t^)h a {\ - (1 - r 2 )x 2 , 1 - (1 - r 2 2 )x 3 , b 2 , b 3 )S t (x 3 )] . (28) 

For the amplitudes of the factorizable diagrams with a c-quark propagator (the first two diagrams in Fig. [3]), we add 
the character "c" in the subscript to distinguish them from those with a light quark propagator. Because of current 
conservation, the factorizablc annihilation diagrams of B —> PP decay mode for (V-A)(V-A) and (V-A)(V+A) 
operators cancel each other. Amplitudes for these diagrams are given as follows: 

• Factorizablc annihilation diagrams for (V-A)(V-A) operator 

Fac L Mt)) - . (29) 

• Factorizablc annihilation diagrams for (V-A)(V+A) operator 

F ac R ( a i(t)) - F a L c L (a t (t))=0. (30) 



Similar to the factorizablc annihilation diagrams, the nonfactorizablc annihilation diagrams are also divided into 
two types (see the last two diagrams of Figs. [2] and [3]), depending on whether the cc are generated from the effective 
weak vertex. Because the c quark in the charmed meson carries most of the energy, these two types of nonfactorizablc 
diagrams are expected to have different scales. Additionally, because the momentum fraction Xi(i = 2,3) is defined 
on the light quark in the charmed mesons, these two types of nonfactorizablc annihilation diagrams also have different 
expressions. The amplitudes of the diagrams with cc pair generated from a hard gluon (the last two diagrams in 
Fig. [2]) are given as 

• Nonfactorizablc annihilation diagrams for (V-A)(V-A) operator 

[2 f 1 f 1/A 

Fani^)) = 167ry -C F m% j [dx] J bidbib 2 db 2 (j} B {xi,bi)4>M 2 {.x 2 )4>M 3 {x?,) 

x (r 2 r 3 (x 2 + x 3 -4)+x 3 - l)E an (t^)ai(t^)h^ (x it h) 
~{r 2 r 3 {x 2 + x 3 - 2) + x 2 - l)E an (t^)ai(t^)h^(^,bi) 

• Nonfactorizable annihilation diagrams for (V-A)(V+A) operator 

[2 f 1 f 1/A 

Fan ^i^)) = 167ry-CFm| J [dx] J bidbib 2 db 2 (j)B{xi,bi)(j)M 2 {x 2 )(l)M 3 {x 3 ) 

x -(r 2 (x 2 + 1 ) - r 3 (x 3 + l))E an {t^)aS { ^)h^iM) 
+ {r 2 (x 2 - 1) - r 3 (x 3 - l))E an (t^)a t (t^)h^(x t ,b z ) . 



(31) 



(32) 



• Nonfactorizablc annihilation diagrams for (S-P)(S+P) operator 

r-l 



2 r L r 1 / A 

F an( a i( t )) = 16tt^ -C F m% J [dx] j bidbib 2 db 2 (j)B{xi,bi)<t>M2{x2)<t>M 3 {x 3 ) 

(r 2 r 3 {x 2 + x 3 -A) + x 2 - l)^^)^^)^^^^) 
r 2 r 3 (x 2 +x 3 -2)+x 3 - l)E an (t^)ai(t^)h^i,bi) 



(33) 



Similar to what we do with the factorizablc annihilation diagrams, the amplitudes with cc pair from the effective 
weak vertex are also distinguished by adding the character "c" in the subscripts. Amplitudes for these diagrams (the 
last two diagrams in Fig. [3]) are given by 



Nonfactorizable annihilation diagrams for (V-A)(V-A) operator 

,i/A 

b 1 db 1 b2db2<pB(xi,bi)4>M 2 (x2)4>M s (x3) 



FL L Mt)) = l<te^C F m% J\dx] 

x [( - r 2 r 3 (x 2 +x 3 + 2)- x 2 )E an (tt ) )a i (tt ) )ht ) M) 
+ {r 2 r 3 (x 2 + x 3 ) + x 3 )E an {t^) ai {t^)h^\xiM) ■ 



(34) 



Nonfactorizable annihilation diagrams for (S-P)(S+P) operator 



16lT\l —CpTllg 



[dx] 



l/A 



bidbib 2 db 2 (f> B {xi,bi)4> M2 {x 2 )(/)M 3 (x 3 ) 



(r 2 r 3 (x 2 +x 3 + 2) + x 3 )E an {t^)a l {t^)h^\x l M) 
((r 2 r 3 + l)x 2 + r 3 x 3 (r 2 + r 3 ))E an (t^) ai (&h^ (x t ,b t ) 



(35) 
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E. Analytic expressions for the decay amplitudes 

There are 10 decay channels for the B — > PP decay mode, which can be divided into two groups: decays with both 
emission and annihilation contributions and pure annihilation type decays. 

• Channels with both emission and annihilation contributions. 

A{B- -+ D°D- s) ) = ^ {v cb V; d{s) [F e LL ( ai ) + F£ £ (d)] + V ub V; d{s) [F a LL ( ai ) + F^{d)\ 

-VtbKn.) ^ LL ^ + °io) + F e L n L (C 3 + d) + F e SP (a 6 + o 8 ) + F™(C S + d) 
+F a LL (a 4 + o 10 ) + F^{d + C 9 ) + F a SP (a 6 + o 8 ) + F^(C 5 + C 7 )}} , (36) 

A(B° -> D+D-) = ^{ Vcb v; d [F e LL (a 1 ) + F^(d)+FtHa 2 ) + F a ^ c (d)} 

-V tb V t * d [F e LL (a 4 + a 1Q ) + F^(C 3 + C 9 ) + F e SP (a 6 + o 8 ) + F e L n R (C 5 + C 7 ) 
+F a L c L {a 3 +a g ) + F^ c (d + do) + F^ R (a 5 + or) + F^ p c (d + d) 
+F a LL (a 3 + a 4 - \a Q - ^a w ) + F^{C 3 + d~ hj Q - ^C 10 ) 

+F a LR (a 5 - \a 7 ) + F^(C 6 - \d) + F a SP (a 6 \a*) + F™(C B - lc 7 )]j , (37) 

A(B° -> D+D-) = ^ {V cb V: s [F e LL ( ai ) + F£(Cx)\ - V tb V t * s [F e LL (a* + a w ) + F e L n L (C 3 + d) 
+F e SP (a 6 + a 8 ) + F e LR (d + C 7 ) + F a Li ( a4 - \a w ) + F^(C 3 - \d) 
+F a SP (a e \a s ) + F a LR (d - \d)]} , (38) 

A(B° S -+ D+D-) = ^ {V cb V; d [F c LL ( ai ) + F e L n L (d)] - V tb V t * d [F e LL (a A + a 10 ) + F e L n L (d + d) 
+F e SP (a 6 + as) + F e L n R (d + d) + F^ L (a 4 - \a 1Q ) + F^{d - \d) 
+F a SP (a 6 \as) + F a LP (d \d)]} , (39) 

A(B° -+ D+D-) = ^ {V cb V; s [F e LL ( ai ) + F e L n L (d) + F a L c L (a 2 ) + F a L n L c (d)} 

~V tb V*{F c LL (ai + aio) + F^{d + d) + F c SP (a 6 + a 8 ) + F e LR (d + d) 
+F a L c L (a 3 + 09) + F^ c (d + do) + Ft R (a 5 + a 7 ) + F a SP c (d + d) 
+Ft L {a 3 + a 4 - ia 9 - ±a w ) + F^(d + d~ \d - \do) 

+F a LR (a 5 - \a 7 ) + F s a p {d - \d) + F a SP (a 6 ~ \as) + F^ R {d ~ \d)]} , (40) 

• Pure annihilation decays. 

A(B° -> D°D°) = ^{V cb V; d [FtH^) + F a L n L c (d)} + V ub V: d [F^ 

-V tb V: d [F^ L {a 3 + ofl) + F^(d + do) + F^ R (a 5 + a 7 ) + F^ p (d + d) 
+F a L c L (a 3 + as) + F^ c (d + do) + F™(a 5 + a 7 ) + F a 5P (d + d)} } , (41) 
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A(B° D+D-) = ^ {V cb V c * d [F a L c L (a 2 ) + F a L ^ c (C 2 )] - V tb V t * d [F^ (a 3 + a 9 ) + F^ C (C 4 + do) 
+i&*(a 8 + a r ) + i^(C 6 + C 8 ) + F a LL (a 3 - ^9) + ^ Q L n L (C 4 - l -C w ) 
+F a LR (a 5 - \a 7 ) + F a s n p (C 6 ^ 8 )| , (42) 

A(B° S ^D°D ) = ^{ Vcb V:AF^ L (a 2 ) + F^ c (C 2 )}+V ub V: s [F a LL (a^^ 

-V tb V; s [F^ L (a 3 + 09) + F a L „ L (C 4 + C 10 ) + F a ifl (a 5 + a 7 ) + F^ (C 6 + C 8 ) 
+^ L (a 3 + 09) + i^ c (C 4 + C10) + Ff c fl (a 5 + 07) + FZ(C 6 + C s )} } , (43) 

A(B° S -+ D+D-) = ^ {V cb V; s [Ft L {a 2 ) + F a L ^ c (C 2 )] - V tb V; s [F^ c L (a 3 + a 9 ) + F Q £ n L c (C 4 + C 10 ) 
+i& Jl (a B + a 7 ) + ^(O, + C 8 ) + ^ Q LL (a 3 - ^9) + F£(C 4 - \c w ) 
+F a LR (a 5 lar) + F*f(C 6 - ic 8 )} . (44) 

There are also 10 decay channels for each category of B —> PV, B —> VP, and B —> VV decays. The decay 
amplitudes of the B —> PV and B — > VP modes can be obtained from the B — > PP decays just by substituting the 
£*( S )/D( S ) meson for the corresponding D* s ^/D* s ^ meson. The factorization formulae for these two decay modes are 
listed in Appendix [B] and [C] respectively. 

The amplitude of B — > VV decay can be decomposed as 

A(e 2 ,e 3 ) = iA N +i(e* 2T -e* 3T )A s + (e^ a0 n»n»e*£ef T )A p , (45) 

where A , including the D wave and part of the S wave component, contains the contribution from the longitudinal 
polarizations A s and A p , corresponding to part of the S wave component and all the P wave component, respectively, 
which represent the transversely polarized contributions, and they have the following relationships with the helicity 
amplitudes (an i in the amplitude is dropped): 

A = A N , A± = A s ± A p . (46) 

For each decay process of B VV, the amplitudes A N , A s , and A p have the same structures as Eq. (|3"6"|) - (|4"4"|) . 
respectively. The factorization formulae for the longitudinal and transverse polarization for the B — > VV decays are 
all listed in Appendix [D] 



III. NUMERICAL ANALYSIS 

The decay widths of B to two charmed mesons decays can be directly derived from the formulas of two-body decays 
in Ref. Q. With the amplitude obtained in Sec. HI the decay widths for the B -> PP, B -> PV, and B -> VP 



decays are given by 



[(l-(r 2 +r 3 )*)(l-(r 2 -r 3 )*)]W 2 _ 

167T771B 
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E i^i 2 - ( 48 ) 



For the B — > VV decays, the decay width is given by 

r = [(l-(r 2 +r 3 ) 2 )(l-(r 2 -r 3 ) 2 )]V 2 

The branching ratio is given by BIZ = Ttb ■ 

The key observables of the decays related in this paper are the CP averaged branching ratios as well as direct CP 
asymmetries (A^p) and mixing induced CP asymmetrics(^4™p x ). Readers arc referred to Ref. [r| for some reviews on 
CP violation. First, we define four amplitudes as follows: 

A f = (f\H\B), A f = (f\H\B), 

A f = (f\U\B), A f = {f\H\B), (49) 
where B meson has a b quark in it and / is the CP conjugate state of /. The direct CP asymmetry Af*p is defined by 

A*> W-I^l 2 (50) 

In neutral B meson decays, if the final states are CP eigen states / = /, the time-dependent CP asymmetry with 
mixing effects present, is defined by 

T(B(t) -+ /) - T(B(t) -+ /) 



A CP (B(t)^f) 



T(B(t)^f)+T(B(t)^f) 

-C f cos(AMi) + A£'p(B -> /) sin(AMt), (51) 



where AM is the mass difference of B meson mass eigenstates. After some calculation, we can get the explicit 
expressions 

\A f \ 2 -\A f \* 



/(mix _ L P ■> 3 ' lrn\ 



\M* + \A f \* 

2Im[lA f A* f ] 
\Af? + \A f 



Since the mixing CP violation in neutral B meson system is negligible in a good approximation, we have 

1 _ „-#M (B) _ V tb V t(d/s) , 

v ~v tb v* d/s) - (b6) 

Our results for CP averaged branching ratios and CP asymmetries are listed in Tables HI UH IIII1 IIV1 and |V] All the 
experimental data are from the Particle Data Group[3] except the ones marked with "BaBar" and "Belle". In Table 
HVl we also list the ratios of the transverse polarizations TZt in the branching ratios for B VV decays, which is 
defined by 



L4 + | 2 + |^L ' 
L4 | 2 + L4 + | 2 + K 



TZ T = 1+1 " r| ^ 1 (54) 



The first errors in our results are estimated from the hadronic parameters: (1) The decay constants of B^ mesons: 
f B = (0.19 ± 0.025)GeV for B mesons and f Bs = (0.23 ± 0.03)GeV for B s meson; (2) The shape parameters in 
Br s \ meson wave functions: uji, = (0.40 ± 0.05)GcV for B meson and uJb B = (0.50 ± 0.05)GcV for B s meson; (3) 
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The decay constants and the shape parameters in the wave functions of charmed mesons, which are given in the 
last paragraph in Sec. Ill CI The second errors are from the not known next-to-leading order QCD corrections with 
respect to a s and nonperturbative power corrections with respect to scales in Eq.(7), characterized by the choice of 
the Aqcd = (0.25 ± 0.05)GeV and the variations of the factorization scales shown in Appendix [A"l The third errors 
are brought in by the CKM matrix elements, which arc given as [2^] 

\V cb \ = 0.041 171858? S. \V cd \= 0-225 OBlglggg g , \V„\= 0-973 4718:8882, 

\v ub \ = o.oos 5i8;888 11, IK,| = 0.974 ut^zil ,\v us \ = 0.225 718:88!!, 
\v tb \ = 0.999 14618:888 8^, W«\ = o.oos 59i8:888S, 1^1 = 0.040 4i±g : gg? ?| , 

7 = (67.8l|-j)°, ^ = (21.5818:^)°. (55) 
The other input parameters are [3| 

G F = 1.16639 x 10- 5 GeV~ 2 , 

t b - = 1.639 x 10~ 12 s/ft . t b» = 1-530 x 10~ 12 s/ft , r B o = 1.478 x 10~ 12 s/ft , 
m B = 5.28GcV , m Bs = 5.366GcV , m D = 1.87GeV , m Da = 1.97GeV , 

m D , = 2.01GcV , m D . = 2.11GeV , h = 6.582119 x 10~ 25 GeV s . (56) 

Because in the direct CP asymmetries the errors arising from the CKM elements are very small, we neglect them. 
In the B — > VV decays, the ratios of the transverse polarizations' contributions (TZ^are not very sensitive to the 
parameters listed above. The next-to-leading order corrections on r occur at the r 2 = 0.13 order and thus the errors 
from the higher orders of r are very small except for TZt- This is confirmed at the numerical calculations. Therefore 
we only keep these errors in TZt and neglect them in other physical quantities. We will talk about the errors of these 
ratios later. 

The first 6 channels in each of Tables Q] HH IIIII IIVI and [V] receive contributions from both emission diagrams and 
annihilation diagrams; while the last 4 channels in each table are pure annihilation processes. In order to make our 
discussions easier, we give a number to each channel in the beginning of each line in the tables. 

Compared with the tree operators, the penguin operators give very small contributions because of the severe sup- 
pression of the Wilson coefficients. By calculating the ratio of the branching fraction with only penguin contributions 
and that with all contributions in the same channel, we estimate how much the penguin operators contribute. Our 
results show that the penguin operators contribute 0.1% — 0.2% in those channels with a W emission contribution, 
and contribute 0.3% — 0.7% in those pure annihilation processes. Thus it is enough to pay our attention only to the 
tree operators in the following for the investigation of the branching ratios. Different from the counting rules of the 
PQCD calculation of B to two light mesons decays, the nonfactorizable emission diagrams may give large contribu- 
tions because the asymmetry of the two quarks in charmed mesons can not make the two diagrams nearly cancel each 
other. However, from Eq. (|36p -Eq.(|44 |) . one can find that the contributions of the nonfactorizable emission diagrams 
are suppressed by the small Wilson coefficient C\. Since the charm quark is heavier than the u, d, s quark, the gluons 
in Fig. [3] are softer than those in Fig. [5J This indicates that the diagrams in Fig. [3] will give larger contributions than 
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TABLE I: CP averaging branching ratios (unit: 10 4 ) and the CP asymmetries for B — > PP decays. 





BTZ 




^CP 


Channels 


Exp. This work 


Exp. This work 


Exp. This work 


1 B- -> D°D" 

2 B~ — > D°D- 

3 B° -> D + D~ 

4 B°-> 

5 B° -» D^-D" 

6 B° a ^D+D7 


4.2 ±0.6 3.9± 2 ; 9 9 ±™.2 
103 ± 17 95±f 6 i«t 2 6 
2.11 ±0.31 S^^oi 

74 ±7 89±S±^±3 

r, o +1.4+0.7+0.1 
z - z -l. 0-0. 7-0.1 

110 ±40 55±H±1S±| 


-13 ±14 ±2 o.6t° :t-°oi 

~ -10 -3 

ll±22±7[Babar] O^t^to! 
-91 ± 23 ± 6 [Belle] 

n r + 0.1 + 0.2 

u - d -o. 0-0.1 


-o.8i ±0.29 -o.73i°;»ro: oii o: o2 


7 B°^D°D° 

8 B° -)■ D+.DJ 

9 B°^D°D 

10 B° -> D + D~ 


<0.6[BaBar] O^St^S 
< 0.36[Belle] 0.35±g||± ;?^;SI 

r n + l. 7+1. 0+0.1 
°- u -l. 5-1. 2-0.3 

r 9+1.5+0.7+0.1 
°- z -l. 9-1.4-0.3 


r q+0.2+0.0+0.2 
J.O_ 2.7-3.3-0.3 

9 q+0.5+0.8 
0.4— 0.4 

n 9+O.I+O.I 
u - z -0. 0-0.0 


n 74+0.00+0.00+0.02 

u - '^-0.01-0.01-0.02 

n 7 o+0. 00+0. 00+0. 02 
u - ' •J-0. 00-0. 01-0. 02 

~ 10~ 3 



TABLE II: CP averaging branching ratios (unit: 10 4 ) and the CP asymmetries for B — > PV decays. 





BTZ 


A 


%>(%) 


Channels 


Exp. 


This work 


Exp. 


This work 


1 


B~ - 


-> D°D*~ 


3.9 ±0.5 


o fi +2. 6+0. 7+0.1 
— 1.7 — 1.0 — 0.2 




n 1+0.4+0.1 
u - A -0. 1-0.1 


2 


B~~ - 


-> D°D* a - 


78 ±16 


QQ+64+20+2 
° 3 -42-24-5 




~ -10~ 3 


3 


B° - 


> D+D*- 


6.1 ± 1.5 


q 9+2.4+0.5+0.1 

■J-^-i. 5-0.8-0.2 


-6 ±9 


~ 10 -2 


4 


B°- 


> d+d;- 


76 ±16 


oq+61+17+2 
°°-39-23-5 






5 


m- 






9 i +1.3+0.7+0.0 
z - l -0. 9-0. 7-0.1 




n i+O.O+O.O 
u - 1 -0.1-0.0 


6 


m- 


* Ds+^r 




/(0+31+15+1 
48_ 21 _ 15 _3 






7 


B°- 




< 2.9[Babar][25j (4.6±^;f ±li±S.l) X 10~ 2 




A 1+1.3+0.0 
l-r-4.4-2.9 


8 


B°- 




< 1.3[BaBar][24j (3.5±J;;f ±l:f ±gi ) x KT 2 




n c+o.i+1.7 

u - d -0.3-0.7 


9 


m- 






fl oq+0. 41+0. 32+0. 01 
u - M -0. 24-0. 19-0. 04 




+0.1+0.1 
u -^-0. 2-0.1 


10 B° - 


> D + D*~ 




n r, A +0.23+0.24+0.01 
'^-0. 29-0. 23-0. 04 







those in Fig. [2] It is confirmed by our numerical results. However, these contributions are still much smaller than 
those from the factorizable emission diagrams. Thus, the branching ratios of the first 6 channels in Tables U [Til Hill 
and IIVI arc dominated by the factorizable emission diagrams. 

Because the factorizable emission diagrams are dominate, the amplitude of the first 6 channels in each table should 
be nearly proportional to the product of a decay constant and a B — > D transition form factor. Based on this physical 
picture, we can have the following simple conclusions: 
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TABLE III: CP averaging branching ratios (unit: 10 4 ) and the CP asymmetries for B — > VP (characterized by B — > V form 
factor) decays. 



BTZ 



,(%) 



Channels 



Exp. 



This work 



Exp. 



This work 



1 


B~ - 


-> D*°D~ 


2 


B~ - 




3 


B°- 


> D* + D- 


4 


B° - 




5 


B° s - 


* d:+d- 


6 


B°- 




7 


B° - 


■4 D*°D° 


<S 


B° - 




9 


m- 


■+ D*°D° 


10 B° - 





6.3±1.4±1.0[BaBar][25] 
84 ± 17 
8.8 ± 1.6 
83 ±11 



4.8 



+3.4+1.1+0.1 



-2.3-1.4-0.3 

11Q+ 94 + 27 + 2 
i±t '-56-34-7 

-I e+3. 5+0. 9+0.1 
y: - u -2. 1-1. 1-0.3 

1 i 9 +86+26+2 
llz -53-32-6 

o 7 +l. 7+0. 9+0.1 
z - ' -1.1-0.9-0.1 

7A+44+19+1 
'"-31-21-4 



-0.5+°- 
10 



0.2-0.3 

3 



-0.6 



+0.0+0.1 
-0.1-0.2 



-0.4 



+0.0+0.1 
0.0-0.1 



< 1.3[BaBar][24] 



0.21 
0.25 



+0.06+0.03+0.00 
0.08-0.05-0.01 



-0.08+0.06+0.01 
-0.08-0.08-0.01 



4.3; 

4.4" 



-1.3+0.8+0.1 
-1.3-1.1-0.2 



-1.4+0.9+0.1 
-1.3-1.2-0.2 



i O+0.3+0.7+0.1 
" J — a -l. 2-1.4-0.1 

~ -10~ 3 



TABLE IV: CP averaging branching ratios for B — } Vl^unit: 10 4 ) and the ratios of the transverse polarizations' contribution. 



BIZ TZt 



Channels 


Exp. 


This work 


Exp. 


This work 


1 B~ - 


-> D*°D*- 


8.1 ± 1.2 ± 1.2[Babar][25] 


f. o+S. 0+1. 5+0.1 
°-°-3.2-2.0-0.4 




0.4571" 


2 B~ - 


■+ D*°D*~ 


175 ± 23 


1 «1 +139+41+3.5 
lol -95-53-10 




0.481°;- 


3 B°- 


> D*+D*- 


8.2 ±0.9 


R n+4. 8+1. 1+0.1 
u - J -3.0-1.6-0.4 


0.43 ±0.08 ±0.02 


0.46^ri4 


4 B°- 


> JJ* 4 -^- 


179 ± 14 


-IOO + 130+39+3.2 
loc5 -88-48-9.6 


0.48 ± 0.05 


0.48^11 


5 B° B - 


+ D*+D*- 




o q+2. 6+1. 2+0.1 
J - M -1.9-1.3-0.2 






6 m- 


+ o: + or 




QQ +72+ 26+1.9 
ay -54-29-5.6 




47+°' 15 

U ' 4 ' -0.15 


7 B°- 




< 0.9[Babar][25] 


n 1 c+0. 05+0. 03+0. 00 
u - io -0. 04-0. 03-0. 01 




n 47+0-35 

u -*'-0.29 


8 B°- 


* £>: + £>r 


< 2.4[Babar][24j 


n 1 Q+0. 10+0. 06+0. 00 
u - ly -0. 07-0. 05-0. 01 




0.571°;!? 


9 m- 






9 o+l.l+0. 7+0.1 
z -°-0.8-0.6-0.2 




0.481^27 


10 B° s - 


+ D *+ D *- 




o i +1.0+0.9+0.1 
J - 1 -0.8-0.7-0.2 




0-49± ;i| 



1. In each table the channels 1 and 3 should have similar branching ratios because they have the same CKM matrix 
elements for the factorizable emission diagrams and similar transition form factors for isospin symmetry. For 
the same reason, channels 2 and 4 should also have similar branching ratios. 

2. The branching ratios of channels 4 and 6 indicate that the B s — > Ds transition has a little smaller form 
factor than B — > D^*' transition. The reason is that the antistrange quark in the B s meson has a little larger 
momentum fraction than the d quark in the B meson due to the SU(3) breaking effect Q. In Q, the 
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TABLE V: CP asymmetry and the ratios of P-wave contributions in branching ratios for B — > VV decays. 



>(.%) 



1 mix 



Channels 



Exp. This work 



Exp. 



This work 



Ri 



1 B- 

2 B~ 

3 B° 

4 5° 

5 B° 

6 5° 



► D*°D*~ 

> D*°D*- 
D* + D- 
D*+DT 
Dl+D- 
Dl+DZ- 



n o+O.O+O.O 
u -^-0. 1-0.1 

~ -10~ 3 
2 ±10 ~ -10~ 2 



-0.01-0.03-0.02 



0.17 
0.16 
0.16 

0.14 
0.17 



7 B° -> D* U D 

8 B 

9 B° B 

10 B° -> D* + D 



o n*o 



-3.4 



+0.4+0.4 



• d: + dt 

D*°D*° 



-0.4 



+0.1+0.3 



-0.73 
-0.68 



0.2 



+0.0+0.0 



+0.03+0.23+0.01 
-0.05-0.14-0.01 

+0.03+0.27+0.01 
-0.06-0.17-0.01 

- 10~ 3 



0.24 
0.32 
0.25 
0.29 



B s — » D s transition is investigated with the light-cone sum rules, and a similar branching ratio for B s — > D~ 
is obtained under the factorization assumption. This means the PQCD and the light-cone sum rules have the 
similar B s — > D s transition form factors. 

3. The first 6 B — > PP decays in Table Q] and the corresponding 6 B — > PV decays in Table [TT| have the same 
transition form factors, respectively, as well as the similar decay constants between D meson and D* meson. 
Thus their branching ratios should also be similar. However, such phenomena arc not expected in B — > VP 
decays and B — > VV decays, because in addition to the longitudinal polarization's contributions, the B — > VV 
decays also receive large contributions from transverse polarizations. 

From Tables HI HU IIII1 and HVl one can find that most experimental branching ratios agree with our conclusions in 



2G 



27| also investigate the decays of B to 



the above paragraphs very well within the errors. The authors in Refs. 
double charmed mesons under factorization assumption, but with different models. Their results also indicate that 
the factorization works well. 

Since the direct CP asymmetry is proportional to the interference between the tree and penguin contributions [lol ] , 
it should be small because of the small penguin contributions as we mentioned above. Our numerical results indeed 
indicate that the direct CP violations are very small. A relatively large direct CP violation appears in the pure 
annihilation decay B° —> D°D° and its corresponding B —> PV, VP, and VV decays. However it is still only several 
percent. Although the experiments give somehow large direct CP asymmetry in some channels, the uncertainty is 
still large. Any large direct CP violation observed in experiments would be treated as a signal of new physics at first. 

The mixing induced CP asymmetry in B decays is almost proportional to the sin 2/3 from Eq.(|53j). if we neglect the 
small contribution from penguin contributions. It should be mentioned that, experimentally the P wave component 
in the amplitudes of B — > VV mode will bring systematic errors in the results of mixing induced CP asymmetry 
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because they will bring a minus sign relative to the S and D wave component. Our results for AJ^p in Table [V] only 
include the S wave and D wave contributions, and in this table we also give the values of R_\_, which is defined by the 
ratio of branching fractions with only P wave component and that with all the contributions. Because the P wave 
contributions are very small in the color allowed tree dominated processes, the experimental measurements are still 
in agreement with our calculations. For the pure annihilation processes, the P-wave contributions arc relatively large 
and therefore these channels may not be good choices for the observation of mixing induced CP asymmetry. 

In Table Hvl we give the ratios of transverse polarizations' contributions in branching ratios. One can find that both 
in the processes with an external W emission and in the pure annihilation decays, the transverse polarizations take 
about 40% — 50% of the contributions, which agree with the present experimental data amazingly well. We should 
point out that these ratios are very sensitive to the terms with power r 2 (r = mu/rnB), although these corrections 
change the other observables only a little. With the r 2 corrections absent, the ratios of transverse polarizations for 
the channels with an external W emission are about 20%, and those for the pure annihilation channels are 0, because 
the transverse polarization's contributions for these channels are at the power of r 2 . For the sensitivity of these ratios 



30(, the authors 



to the power correction terms, we vary the variable r by 20% for an error estimation in Table IIVI In 
obtain the values for the ratios ~ 50% for the external W emission processes simply by means of kinematics under 
the naive factorization, which agree with our results. For the pure annihilation decays, the transverse polarizations 



30( think the ratios of the transverse polarizations are 



are suppressed by r 2 , which is the reason why the authors in 
very small. However, our calculation show that, with the r 2 terms included, these ratios increase to about 50%. This 
means that the polarization fractions are quite sensitive to the power corrections although they are not sensitive to 
the higher order QCD corrections etc. The future experiments will tell us more about the polarizations in the pure 
annihilation processes. 



IV. SUMMARY 



Although the D meson mass is not very small compared with the B meson mass, factorization can still work in 
the leading order of mo /tub and Aqcd/"1d expansion. Since the PQCD approach can eliminate the end-point 
singularity in the perturbative calculation, we investigate the decays of B to double charmed mesons systematically. 
Both pseudoscalar and vector charmed mesons are included in the final states. We find that the factorizable emission 
diagrams are dominant in the branching ratios. Most of our branching ratios agree with the experimental data, 
which means the factorization assumption works well. However, experimental data show that there are still some 
discrepancies, which means more work is needed both at the theoretical side and the experimental side. 

Our results indicate that the direct CP asymmetries in these channels are very small. Thus, it will be a signal 
of new physics if a large direct CP asymmetry appears. In the decays of B to double vector charmed mesons, the 
transverse polarizations contribute 40% — 50% both in the external W emission processes and in the pure annihilation 
decays, which agree with the present experimental data. We should mention that the correction terms at the power of 
r 2 play an important role in transverse polarizations, without which the ratios for the external W emission processes 
decrease to about 20% and for the pure annihilation decays the ratios are because of the r 2 suppression. 
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Appendix A: scales and functions for the hard kernel 

The variables that are used to determine the scales and the expressions of the hard kernels are defined by 

Pen = m\xYX 2 {\ - r|), 
P e W = m\x 2 {x x {l-r\)-x 3 {\-r\-r\)\ 

Pen = -1Ti B [x 2 (x 3 - + T% (xi{x2 - 1) + X 2 (x 3 - 1) - X 3 ) - X 2 (xi + X 3 - 1)], 

Pan = m%(l - (1 - rl)x 3 -x 2 (l- x 3 (l - r\) - (1 - x 3 )r\)\ 
P£) = m 2 B (l + Xl x 2 (l-r 2 3 )-(l-rl-rl)x 2 x 3 ), 

P$ = m%(-x 3 ti + xi((ti - l)x 2 + 1) + x 3 + x 2 ((x 3 - l)r§ + (r£ - l)x 3 + 1) - 1), 

= m%(l - r\ - r 2 )x 2 x 3 , 
Pin ] = ™lM(r 2 - l)x 2 + 1) - (r 2 2 - l)x 3 + x 2 ((x 3 - l)r 2 3 + (r 2 2 - l)x 3 + 1)], 

P^l c) = <x 2 ((rj + r 2 - l)x 3 (r§ - l)x x ). (Al) 
The scales are determined as 

= max{y / a;2(l - r|)m B / crr , 1/bi, I/&2}, 
4 2) = max{y / xi(l - r|)m B / err , I/61, l/6 2 }, 
t& = max{yjp^|/ err ,yiP^|/ err , I/61, 1/63}, 

45 = max{V^|/ e rr,vi^|/ e rr, l/6l,l/68}. 

4 X) = max{yl - (1 - rl)x 2 m B f m , l/6 2 , 1/63}, 
4 2) = max{y / l - (1 - rj)x 3 m B f c „, V&2, l/&3>, 

i« - max{v^^/err,vl^l/err, l/6l, I/62}, 
4n = max {v^Q/crr, yl^l/err, 1/&1, l/fc}, 

4 lc) = max^ (1 - r\ - r%)x 3 m B /err, I/62, 1/63}, 
4 2c) = max{^/ (1 - ~r\ ~ r|)x 2 m B / crr , l/6 2 , 1/63}, 

4^ = max{V^|/«r, yj |Pin C) |/err, l/6 lf l/6 2 }, 

= maxi-y/^/err, ^ SS/err, 1/&1, l/& 2 }, (A2) 

with / elT varies from 0.75 to 1.25 for an error estimation. 
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The functions of the hard kernels that appear in the factorization formulae are given by 
h e (x 1 ,x 2 ,b 1 ,b 2 ) = K {^xix 2 m B bi) 

x [0(&i - b 2 )K (V^2 m Bh) h (V^Wfly 



h a (x 2 ,x 3 ,b 2 ,b 3 ) 



+9(b 2 - bi)K (y/x^m B b 2 ) Iq (y/x^m B bi)] , 

2 J 



I -H'o 1) Ux 2 x 3 m B b 2 ) 



6{b 2 - b 3 )H^ (y^m B b 2 ) J {y/x^m B b 3 ) 
+9{b 3 - b 2 )H^ (y/x~3~m B b3) J (y/x^m B b 2 ) , 
9(h - b 3 )K (v^fci) /o (V^&s) + #(&3 - 61) 



K {J\P^\b 3 ) forP^>0 



(.;) 



fH^(^\P^'\b 3 ) forP^<0 

- & 2 )i4 1} (v^&i) Jb ( v^fe) + ^2 - 
K (J\pW\bi) for Pi£ > 



1 



Panbl J( 



Panbl 



and the functions that consist of coupling constant and Sudakov factors are given by 

E e {t) = a s (t)exp[-S B (t)-S M2 (t)}, 
E a (t) = a a (t)exp[-S Ma {t)-SAt s (t)], 
E en (t) = a s (t)exp[-S B (t) - S Ah - S M3 \b 2 =bi] , 
E an (t) = a s (t) cxp[- S B (t) - Sm 2 - S M3 \b 3 =b 2 } , 



where 



S B (t) = Sm 2 = Sm 3 = s ( xi^r,bi ) + I 

V2 / 3y 1/6l /x 



^7«(a s (/i)), 



with the quark anomalous dimension 7 g = —a s /ir. The explicit form for the function s(Q, b) is 



s(Q,6) 



2/3i \b 



L ' (- l\ A(2) 



1 



+ 



ln(2g) + 1 ln(26) + 1 



8/3? 



ln 2 (2q) -ln 2 (2&) 



In 



where the variables are defined by 



and the coefficients AW and are 



AC 1 



q = \n[Q/(V2A)}, S = In[l/(6A)], 



_ 33 - 2n f 153 - 19n, 

Pi — — — — j P2 — 



12 



24 



4 f2 v 67 7T 2 10 8 _ , ,1 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



rif is the number of the quark flavors and je is the Euler constant. We will use the one-loop running coupling 
constant, i.e. we pick up the four terms in the hrst line of the expression for the function s(Q, b). 
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Appendix B: factorization formulae for B — > PV(M2 is a pseudoscalar meson and A/3 is a vector meson) 



l/A 



F e (oi(t)) = 8irC F f M3 m B / cteid:E 2 / bidbib 2 db 2 (/) B {xi,bi)<pM 2 {x 2 ) 

Jo Jo 

x -((2a; 2 - l)ra -x 2 - l)E e {t^)a l {t^)h e {x ll x 2 {\ - r 2 )M,b 2 )St{xi) 

(x 2 ,xi(l - rl),b 2 ,bi)St(xi)\ , 



F e Mt)) = , 



.i/A 



F er^( a i( t )) = ^\J^C F m% J [dx] J b 1 db 1 b 3 db 3 (f) B (x 1 ,bi)(f)M 2 {x 2 )(l)M 3 (x3) 
x ( -x 2 r 2 +x 3 )£; 6 (t«)a l (4 1 „))/ l W(x l ,6 I ) 
+ (x 2 r 2 - x 2 + x 3 - l)E en {tf}) ai {tf2)hf2{xiM)\ , 

[2 f 1 f 1/A 

Fen Xa-iit)) = l§ir\j-C F m B J [dx] J bidbib 3 db 3 <j) B (xi,bi)(t) M2 (x 2 )4>M 3 (x 3 ) 



xr 3 



{x 3 - r 2 (x 2 - x 3 ))E en {t^)ai{t^)h^{ Xi M) 



(x 3 + r 2 (x 2 + x 3 ))E en (tW) ai (tW)hW( Xh bi ) 
• l/A 



F T / L (a i (t)) = 8irC F f B m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 4> M2 (x 2 )4> M3 {x 3 ) 

Jo Jo 

x (x 2 - l)E a (tV>)ai(tW)h a (l - (1 - r|)z 3 , 1 - (1 - r 2 3 )x 2 ,b 3 , b 2 )S t {x 2 ) 

+ ( - 2r 2 r 3 x 3 - (x 3 - l))E a {t^) ai (t^)h a (l - (1 - r|)z 2 , 1 - (1 - r§)x 3) 6 2j b 3 )S t {x 3 ) 



-l/A 



F^ p (ai(t)) = 16irC F f B m' B I dx 2 dx 3 / b 2 db 2 b 3 db 3 4> M2 (x 2 )(l)M 3 (x 3 ) 

Jo Jo 

x ra(l - x-,)^^ 1 ))^ (4^)^(1 - (1 - rl)x 3 ,l- (1 - rf)z 2 , b 3 , & 2 )S t (z 2 ) 

+ (2r 2 + r 3 (z 3 - l))K(t£ 2) )ai(f£ 2) )Ml - (1 - r|)z 2 , 1 - (1 - r 2 2 )x 3l b 2 , b 3 )S t (x 3 ) 

,1 ,1/A 

F^ c L (a t (t)) = &ttC f f B m B / dx 2 dx 3 \ b 2 db 2 b 3 db 3 4> M2 {x 2 )cj)M 3 (x 3 ) 

Jo Jo 

x "(r 2 r 3 (l - 2x3) - a:3)^a(ti lc) )ai(*i lc) )ha(a;2, x a (l - r\ - r\), b 2 , b 3 )S t (x 3 ) 
+ (r 2 r 3 + x 2 )E a (t^) ai (t^)h a (x 3 ,x 2 (l - r\ - r 2 3 ),b 3 ,b 2 )S t (x 2 ) , 
F a L c R ( ai (t)) = -F a L c L ( ai (t)), 

[2 f 1 f 1/A 

Fan (°i(*)) = M>ny-C F m B y [dx] y bidbib 2 db 2 4> B {xi,bi)4>M 2 (x 2 )4>M 3 {x 3 ) 

x (r 2 r 3 (a;3 - x 2 ) + x 3 - l^^nH^in^n^i, 6<) 
+ (r 2 r 3 (x 3 - x 2 ) -x 2 + l)E an (tW)ai(tW)h a 2 2(xi,bi) 



(Bl) 
(B2) 



(B3) 



(B4) 



(B5) 
(B6) 



(B7) 



(B8) 
(B9) 



(BIO) 
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[2 f 1 f 1/A 

16n^J-C F m B J [dx] J b 1 db 1 b 2 db 2 (j) B (x 1 ,b 1 )(j)M2{x2)4>M 3 (x 3 ) 

x -(r 2 (x 2 + 1) - r 3 (x 3 + l))E an {t^) ai {t^)h^{x h h) 
+(r 2 (x 2 - 1) - r 3 (x 3 - l))E an (tW)ai{tW)hW(xi, h) , 



(Bll) 



F anc( a i(t)) 



\l —Cpm 



1 ,1/A 

[dx] / bidbib 2 db 2 (f>B{xi,bi)(j>M 2 {x2)4>M 3 (x 3 ) 



Jo 

~{r 2 r 3 (x 2 - x 3 ) + x 2 - l)E an {tV) ai {t^)h^{x h h) 
(r 2 r 3 (x 2 - x 3 ) - x 3 + l)E an (t^)ai(t^)h^(xiM) , 



[2 f 1 f 1/A 

IQttJ -C F m% / [dx] / bidbib 2 db 2 4>B{xi,bi)4>M 2 {x 2 )(f>M 3 (x 3 ) 

x \r 2 r 3 (x 3 - x 2 ) - x 2 )E an {t^)aMan ] )h { ^H^ h) 
+ (r 2 r 3 (x 3 - x 2 ) + x 3 )E an (t^) ai (t^)htH^A) , 



(B12) 



(B13) 



1/A 



I&tt \j^C F m% J [dx] J bidbib 2 db 2 4> B {xi,bi)(t)M 2 {x 2 )4>M 3 {x 3 ) 

x \r 2 r 3 {x 3 - x 2 ) + x^EUC^aMn^H^M) 

~((r 2 r 3 + l)s 2 - r 3 x 3 (r 2 + r 3 ))E an {t^)a^)h^ (x l: h)] , 



(B14) 



Appendix C: factorization formulae for B — s> VP(A'h is a vector meson and M3 is a pseudoscalar meson) 



prion®) 



/•1 rVA 
8nC F f M3 m B I dxidx 2 / hdbib 2 db 2 (j) B (xi,bi)(j)M 2 (x 2 ) 
Jo Jo 

x -((2^2 - l)ra -a?a - l)S e (tW)a i (tW)/i e (xi ) x 2 (l - rf), 61, 6 2 )S t (a: 2 ) 
+r 2 (l + r 2 ) J B e (4 2 ))a l (4 2 ))/ le (x 2 ,x 1 (l-r2),6 2 ,6 1 )5 t (x 1 )l , 



/■l ,1/A 

F e p (a,i(t)) = 16irC F fM a m B / dx x dx 2 / b\db\b 2 db 2 § B {x\, &i)0m 2 (^2) 

Jo Jo 

xr 3 [(r 2 x 2 - l)^ e (4 1 ))a i (4 1 ))/i e (a; 1 , : E2(l - r 2 ), 61, ba)S t (x 2 ) 
-r 2 E e (4 2 ) )oi (4 2) )h e (x 2 , X! (1 - r 2 ) , 6 2 , 6i )5 t (X! )] , 

/2 r 1 /" 1/A 

167ry -C F m B J [dx] J b 1 db 1 b 3 db 3 (j) B {xi,bi)(t)M 2 {x 2 )4>M 3 {x 3 ) 

-( - x 2 r 2 - x a )E b {t^)o i (t^)h^{x i ,b i ) 
x 2 r 2 -x 2 + x 3 - l)E en (tW)ai(t£!)hW(xi, b t 



F e L n L ( ai (t)) 



(CI) 



(C2) 



(C3) 



[2 f 1 f 1/A 

16n^J-C F m B J [dx] J bidbib 3 db 3 (f)B(xi,bi)(l>M 2 (x 2 )<t)M 3 (x 3 ) 

r 2 (x 2 - x 3 ) + x 3 )E m (1&)ai{tV>)hM(xi,bi) 



xr 3 



'2 



(s a + x 3 ~2)-x 3 + 2)E m (tW)a i (tW)hP>(x i , b z 



/•l /-I/A 

&TrC F f B m 4 B / dx 2 dx 3 / b 2 db 2 b 3 db 3 (f>M 2 {x 2 )(pM 3 (x 3 ) 
Jo Jo 

x [(2r 2 r 3 x 2 + x 2 - l) 

x^^ln^tl - (1 - ri)i 3 , 1 - (1 - r\)x 2 , 63, 5 2 )S t (x 2 ) 
+ (l-x 3 )^ (^)a i (^ 2 ))ft a (l - (l-r 2 3 )x 2 ,l- (l-r 2 2 )x 3 ,b 2 ,b 3 )S t (x 3 ) 
= -F a LL (a t (t)), 



p 



r 1 ri/A 

(ai(t)) = l6irC F fBm B / dx 2 dx 3 / b 2 db 2 b 3 db 3 4>M 2 {x 2 )4>M 3 {x 3 ) 
Jo Jo 

x -(2r 3 +r 2 (x 2 - l)) J B a (^ 1 ))a i (4 1 ))/i a (l - (1 - r 2 )x 3 ,l - (1 - r 2 3 )x 2 , b 3 , b 2 )S t (x 
+r 3 (x 3 - l)E a {tW) ai (tW)h a (l - (1 - r§)a;2, 1 - (1 - r 2 2 )x 3 ,b 2 ,b 3 )S t (x 3 )] , 



Co (*(*)) 



Fl P ( ai (t)) 



,1 ,1/A 

8irC F f B m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 (f)M 2 {x 2 )4>M 3 (x 3 ) 
Jo Jo 

x [( - r 2 r 3 - x 3 ) J B a (4 lc ))a i (4 lc ))/i (x 2 , a; 3(l - r 2 - r 2 ), 6 2 , 5 3 )S t (x 3 ) 
-(r 2 r 3 (l - 2.t 2 ) - x 2 )a.j(4 2c) )/ia(a;3, x 2 (l - rf - r 3 ), 63, b 2 )S t {x 2 ) , 

/2 r 1 /" 1/A 

= 1671-y -C F m% J [dx] J b 1 dbib 2 db 2 (f> B (xi,bi)4>M 2 (x2)(t>M3(x3) 
x -{r 2 r 3 (x 2 ~ x 3 ) ~x 3 + l)E an (t^)ai{t^)h^{xi,bi) 
+ {r 2 r 3 (x 3 - x 2 ) -x 2 + l)EM(tW)ai(tW)hV>(xi, &,:)] , 

[2 f 1 f 1/A 

= l&7ry-C F m B J [dx] J b 1 db 1 b 2 db 2 (t)B(xi,bi)^M 2 (x 2 )(t>M 3 (x 3 ) 

x (r 2 (x 2 + 1) - r 3 (x 3 + l))E an (tW)ai(tU)h£>(xi, h) 
-(r 2 (x 2 - 1) - r 3 (x 3 - l))E an (t^2)ai(t^)h^M) , 

[2 f 1 f 1/A 

= 167i-y -C F m% J [dx] J bidbib 2 db 2 <j) B {xi,b 1 )(j)M 2 {x 2 )4iM 3 {x 3 ) 

x (r 2 r 3 (x 3 - x 2 ) -x 2 + \)E an (t { V) ai {t^)h^{x i: h) 
~{r 2 r 3 (x 2 - x 3 ) -x 3 + l)E an {t^) ai {t^)h^{xiM)\ , 
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FanMt)) 



16n\/ — C F m 

->o Jo 



l/A 



Ff n P Mt)) = 



[dx] / bidbib 2 db 2 (f>B{xi,bi)4>M2(x2)4>M 3 {x3) 
Jo 

r 2 r 3 (x 3 - x 2 ) - x^E^t^a^an^Kx^ h) 
-(r 2 r 3 (x 2 - x 3 ) - x 3 )E an (t^)a t (t^)h^\ Xl ,b t ) , 

[2 f 1 f 1/A 

167ry -C F m% J [dx] J b 1 db 1 b 2 db 2 ^B(xi,b 1 )(j) M2 (x2)^M 3 (x 3 ) 

x -{r 2 r 3 (x 2 - 13) - x 3 )E an {t^)ai(t^ ) )h^ ) (x i ,bi) 
{{-r 2 r 3 - l)x 2 + r 2 r 3 x 3 )E an (tt ) )a i (&h^(x i ,b i ) 



(C13) 



(C14) 



Appendix D: factorization formulae for B — > VV 



1. Longitudinal polarization 



.l/A 



F e LL (a i (t)) = 8TrC F f M3 m% / dx x dx 2 \ hdbib 2 db 2 (/) B (x 1 ,b 1 )<l>M 2 (x 2 ) 

Jo Jo 

x (-l-x 2 +r 2 (2.T 2 -l))^ e (4 1 ))a i (4 1 ))/i e (xi,x 2 (l-r2),6i,6 2 )^(a ;2 ) 
~r 2 {l + r 2 )E e {t^)a l {t^)h e {x 2 ,x 1 {l-rl)MM)S t {x 1 )\ , 
F e SP Mt)) = 0, 

[2 f 1 f 1/A 

F en = 16ny-C F m% y [dx] y b 1 db 1 b 3 db 3 (j) B {xi,b 1 )(j)M2{x 2 )(j)M 3 {x 3 ) 

x (-x 2 r 2 - X3)£b(4nVi(4M^A) 

~{x 2 r 2 ~x 2 + x 3 - l)E m (tW)ai(tW)hW(Xi,bi) , 

/2 f 1 f 1/A 

Fen Xa-iit)) = 167ry -C F m% J [dx] J bidbib 3 db 3 <f) B (xi,bi)(t) M2 (x 2 )4>M 3 (x 3 ) 

xr 3 -(r 2 (x 3 +x 2 ) ~ x 3 )E en {t^)a l {t^)h { f ^{x l ,b i ) 
-{x 3 + r 2 (x 2 - x 3 ))E en {tM)ai{t¥2)hW( Xubi ) , 

.1 .l/A 

Fa L (a,i(t)) = %irC F f B m B I dx 2 dx 3 \ b 2 db 2 b 3 db 3 4>M 2 {x 2 )(j>M 3 {x 3 ) 

Jo Jo 

x (-x 2 + l)E a {t^) ai {t^)h a {l - (1 - r|)x 3 , 1 - (1 - r\)x 2 , b 3 , b 2 )S t {x 2 ) 
-(1 - x 3 )E a {t^) ai {t^)h a {l - (1 - r2)x 2 , 1 - (1 - r^a*, 6 2l 63)^(3:3)] , 
F a LR (a t (t)) = F a LL ( ai (t)), 

,1 .l/A 

F a( a i{t)) = 167rCir/ B TOB / cte 2 cfe 3 / b 2 db 2 b 3 db 3 (j)M 2 {x 2 )(j>M 3 {,x 3 ) 

Jo Jo 

x r 2 (x 2 - l) J E a (4 1 >)a i (iW)/i a (l - (l-ri)a:3,l- (1 - r\)x 2 , b 3 , b 2 )S t (x 2 ) 
+r 3 (x 3 - l)E a (tW) ai (tW)h a (l - (1 - r\)x 2l 1 - (1 - r 2 2 )x 3 , 6 2 , b 3 )S t {x 3 ) 



(Dl) 
(D2) 



(D3) 



(D4) 



(D5) 
(D6) 



(D7) 



FL L Mt)) 



F a S n P J ai (t)) 



0. 

^"(oi(*)) = o, 



[2 f 1 f 1/A 

^^-C F m%J [dx] J hdbib 2 db2<t> B (xi,b-i)(l>M 2 {x2)<t>M 3 {x3) 

( - r 3 r 2 (x 2 + x 3 ) ~x 3 + l)£ on (tW)o i (*^)/iW(i <> 6 i ) 
- r 2 r 3 (x 2 + x 3 - 2)~x 2 + l)E m {tW)ai(tW)hW(xi, h 



1/A 



M>7ry^C F m B J [dx] J bidbib 2 db 2 <f) B (xi,bi)(l>M 2 (x 2 )<f)M 3 (x 3 ) 

x -(-r 3 (l + 13) + r 2 (l + x 2 ))S a „(tW)a,(tW)/iW(x 4 , & 4 ) 
+(r 2 (z 2 - 1) - r 3 (.T 3 - l))E an (t^)ai(t^)h^(xi,bi 



1/A 



16Tr\l-C F m B [dx] bidbib 2 db 2 (j> B {xi,bi)4>M 2 {x 2 )(j)M 3 (x 3 ) 
Jo Jo 

- r 2 r 3 (x 2 + x 3 ) - x 2 + l)E an {tW)ai(t£2)hW(xiM) 
- r 2 r 3 (x 2 + x 3 - 2) - (x 3 - l))E an (t^)a l (t^)h^(x i ,b i )} , 



16tt\I — C F m B I [dx] 



1/A 



bidbib 2 db 2 (f) B {xi,b 1 )4> M2 (x 2 )4>M 3 (x 3 ) 



( - r 2 r 3 (x 2 +xs-2)- x 2 )E an (t^) ai (t^)h^ (x^) 
+ ( - r 2 r 3 (x 2 + x 3 ) - x 3 )E an {t^) ai {t^)h^\x u &;)] , 

[2 f 1 f 1/A 

lQny -C F m B J [dx] J bidbib 2 db 2 <i> B {xi,b 1 )(j)M 2 {x 2 )4iM 3 {x 3 ) 

x - ( - r 2 r 3 (x 2 +x 3 -2)- x 3 )E an (t<$ K (t^)h^ (x s , &,) 
+ ((-r 2 r 3 - l)x 2 - r 3 r 22 ;3)i? a „(e c) )«»(C c) )^ c) (^, h)} , 



2. Transverse Polarization 



1/A 

F^ L ' s {ai{t)) = 8TrC F f M3 m B I dxidx 2 I bidbib 2 db 2 (j) B (xi,bi)(t)M 2 (x 2 ) 



xr 3 



(r 2 (x 2 + 2) + l)E^) ai {t^)h e {x u x 2 (l - r 2 3 )MM)St{x 2 ) 



-r 2 E e (t^) ai (t™)h e {x 2 ,x 1 (l-rl),b 2 ,b 1 )S t (x 1 ) 



F^ L ' p (ai(t)) = 8nC F f M3 Tn% I dxidx 2 / b 1 dbib 2 db 2 (j) B {x 1 ,bi)4'M2{x 2 ) 



1/A 



xr 3 



(r 2 x 2 -l)S e (4 1 ))a 4 (4 1) )^(xi,a;2(l-r2),6 1 ,6 2 )5 t (x 2 ) 
~r 2 E e (4 2 ) )cn (4 2) ) h e (x 2 , x x (1 - r\ ) , b 2 , h )S t {x x )] , 
if p ' 8 (a,(i)) = F e sp *(a i (t))=0 J 



[2 f 1 f 1/A 

5 (a,;(t)) = 16n\J-C F m B J [dx] J h dbi b 3 db 3 (t>B(xi, 61)^2(^2)^3(^3) 

-XsEbit^aM^h^ix^h) + (1 + x 3 + r 2 (2x 2 - 2x 3 + 1)) E en (t^) ai (t^)h^ ( 



xr 3 



1/A 



Fen' P ( a i(' t )) = 16^^/1^771^ y [cfe] y b\db\bzdbz(j>B{xi, bl)4>M 2 {X2)4>M 3 (X3) 

xr 3 [-^(^K^M^i, 6,) + (1 + x 3 - r 2 )S e „(t( 2 „))a l (4 2 n ))/ 1 ( 2 „)(x J , 6,-) 



[2 f 1 [ 1/A 

l6-Ky-C F m%J [dx] J hdbibsdb 3 (f>B(xi,bi)4>M 2 (x2)(f>M 3 (x3) 

x -(x 2 r 2 - a; 2 r 2 + r^ 3 )S en (tW)o i (*g l ))/»W(i <> 6 j ) 
-(.x 2 r 2 - x 2 r 2 - 

/2 r 1 r 1/A 

l&K\J-C F m B J [dx] J b 1 dbib 3 db 3 <t>B{xi,bi)(t>M2(x2)<t > M s (x3) 

x -{x 2 r\ - x 2 r 2 - rlx 3 )E en {t^2)ai{t^2)h { ^,{xi, h) 
-{x 2 r\ - x 2 r 2 +rlx 3 )E en (t^)a i (t^)h { ^{x i ,bi 



/•l /-1/A 

8>'KC F f B m B / dx 2 dx 3 I 1 
^0 Jo 



xr 2 r 3 



b 2 db 2 b 3 db 3 (f)M 2 (x 2 )4>m 3 (x 3 ) 



(2 - ^^(tWH^Ml - (1 - r£)x 3 , 1 - (1 - r 2 )a; 2 ,63,6 2 )5 t (x 2 ) 



+ (.t 3 - 2)£ B (f< a >)ai(t( 3 >)/i (l - (1 - r 2 )x 2 ,l - (1 - r 2 )^, 6 2 , b 3 )S t (x 3 ) 

/■l fl/A 

L ' P ( a i(*)) = 8nC F f B m 4 B / dx 2 dx 3 / b 2 db 2 b 3 db 3 4>M 2 (x 2 )4>M 3 {x 3 ) 

Jo Jo 

xr 2 r 3 [x 2 E a (tW)ai(tP)h a (l - (1 - r 2 )x 3 , 1 - (1 - r 2 )a; 2 , 63, & 2 )5 t (x 2 ) 
+.T 3 £ a (4 2) H(4 2) )Ml - (1 - r 2 )^2,l - (1-^3,62,63)^3) , 

^(oi(*)) = -F a LL *(a t (t)), 

r-l /.1/A 

Ff^Mt)) = l6TTC F f B m 4 B / da; 2 da;3 / b 2 db 2 b 3 db 3 (j) M2 (x 2 )(j)M 3 (x 3 ) 

Jo Jo 

x -^(^'M^Ml - (1 - r*>3, 1 - (1 - r 2 )x 2 , 63,62)^2) 
-r 2J B (4 2) )ai(4 2) )^(l - (1 -4)x 2 ,l - {1 - r 2 2 )x 3 ,b 2 ,b 3 )St(x 3 )\ , 

,1 j-1/A 

^a L c L ' s ( a i(0) = 8wC F fBm% / dx 2 dx 3 / b 2 db 2 b 3 db 3 4>M 2 (x 2 )4>M 3 (x 3 ) 

Jo Jo 

x -r 2 (r 2 - r 3 (x 3 + l^Ct^M*^)/^, x 3 (l - r\ - r 2 ), b 2 , b 3 )S t (x 3 ) 
+r 3 {r 3 - r 2 {x 2 + l))E a (t^)a t (t^)h a (x 3 ,x 2 (l - v\ - r 2 3 ),b 3 ,b 2 )S t (x 2 )] , 



i& L,p M*)) 



&7rC F fBm B / dx 2 dx 3 / b 2 db 2 b 3 db 3 (j)M 2 (x 2 )(l>M 3 (x 3 ) 
Jo Jo 

x ra(ra - r 3 (z 3 - l))£ a (4 lc) W(4 lc) K(^2, z 3 (l - r| - r 3 2 ), 6 2 , b 3 )5 t (a; 3 ) 
+r 3 (r 3 -r 2 (z 2 - l))K(t£ ac) )oi(t<^)/i„(a:3, z 2 (l - r 2 2 - r 3 2 ), b 3 , b 2 )S t (x 2 ) 



?LL,. 



M*)) 



FL R ' P Mt)) = F^{ ai {t)) 



16n\ —C F m, 



[dx] 



1/A 



bi dbi b 2 db 2 4>b {x\ , b\ )4>m 2 (x 2 ) 4>m 3 (2:3 ) 



(z.r 2 - 2r 2 r 3 + rtx 3 )E an {t^)ai{t^)h^{xi, h) 
((x 2 - l)r 2 + r 2 3 (x 3 - l))E an (t^)ai(t^)h^uh) 



[2 f 1 f 1/A 

167ry -C F m% J [dx] J hdbib 2 db 2 ^ B (xi,b 1 )(f>M 2 (x 2 )(t)M 3 (x3) 

x (r 2 3 x 3 - r\x 2 ) E an (i£> )a, t (x, t , h) 

+ (r 2 (x 2 - 1) - r 2 (x 3 - l^an^M^lSfo. > 

/2 r 1 /" 1/A 

lQ-K\j-C F m% J [dx] J b 1 db 1 b 2 db 2 4> B {xi,b 1 )(f)M 2 (x2)4>M 3 (x 3 ) 

x (ra (x 2 + 1) - r 3 (a* + l))£ a „ (t£>)a< (i$ , 6,) 
-(r 2 (aj a - 1) - r 3 (x 3 - l))E an (t^)ai(t^)h^(xi, Jh)] , 



i^ p (ai(t)) = 



/2 A 1 /" 1/A 

167ry -C F m% J [dx] J b 1 db 1 b 2 db 2 (f) B (x 1 ,b 1 )(t) M2 (x2)(f>M 3 {x 3 ) 

x - {x 2 v\ - 2r 2 r 3 + r 2 x 3 )E an (i« )a, (i « (*< , 6,) 
+ ((x 2 - l)r 2 2 + r 2 (x 3 - l))^^))^^))^)^,^)] , 

/2 r 1 r 1/A 

-C F m B J [dx] J b 1 dbib 2 db 2 4> B {xi,b 1 )4>M 2 {x 2 )4>M 3 {x 3 ) 
x (r 2 x 2 - rlx 3 )E an (t^2)aM2)h^i,h) 
-{r\{x 2 - 1) - r 2 3 {x 3 - l^Eanit^aiit^h^ixiA)} , 

[2 f 1 f 1/A 

167ry -C F m B J [dx] J hdb 1 b 2 db 2 (t) B {x 1 ,bi)(j)M 2 {x2)(t>M 3 (x 3 ) 

\{x 2 - l)r 2 + 2r 2 r 3 + r\{x 3 - l))E an {t^) ai {t^)h^ {x u h) 
x 2 r\ + x^EU&a^&h^Kx^bM , 
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?SP,p 



l/A 



I / bidbib 2 db24>B(xi,bi)4>M 2 (x2)4>M 3 (x 3 ) 
o 



x '-(rl(x 2 - 1) - r 2 (.* 3 - l))^^)^^)^^,^) 



.l/A 



l&TT\j^C F m% J [dx] J bidbib 2 db 2 4>B(xi,bi)(l)M 2 (x2)4'M 3 {x3) 
x \{x 2 - l)r| + 2r 2 r 3 + r 2 (x 3 - l))E an {t^) ai {t^)h^\x^ h) 
+ (r|(x 2 - 1) - rl(x 3 - l))E an {t^) ai {t^)h^\x % M)] , 



M*)) = 



l/A 



I / bidbib 2 db 2 (j)B(xi,bi)4>M 2 (x 2 )4>M 3 (x3) 
o 



x /r 2 2 (x 2 - 1) - r 3 2 (a;3 - l))E an {t^y)a^)h^ {x u h) 
-(r 2 x 2 -r 3 2 x 3 )ii; a „(e ) )^(e c) )^ c) (^,^)" , 
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